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Aiming at the link between confinement and chiral symmetry the Polyakov loop represented as a 
spectral sum of eigenvalues of the Dirac operator was subject of recent studies. We analyze the 
volume dependence as well as the continuum behavior of this quantity for quenched QCD using 
staggered fermions. Furthermore, we present first results using dynamical configurations. 
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1. Motivation 



Lattice simulations suggest that the chiral phase transition and deconfinement phase transition 
appear at the same temperature. It is believed that there is a connection between both phase tran- 
sitions. While for the chiral phase transition we have a well established picture of the symmetry 
breaking mechanism the picture of the deconfinement phase transition remains unclear. Although 
there has been much progress in the last years the final link connecting both phase transitions is 
still missing. 

The order parameter of the chiral phase transition is the chiral condensate {qq). Banks and 
Casher [p related the chiral condensate to the eigenvalue density p of the Dirac operator near zero, 

{qq) = -7Zp{X=Q). (1.1) 

Recently, Gattringer established a formula which relates the eigenvalues of the Dirac operator to 
the Polyakov loop P, the order parameter of the deconfinement phase transition in the quenched ap- 
proximation. This relation provides a natural link between the chiral condensate and the Polyakov 
loop via the eigenvalues of the Dirac operator. The hope is to obtain some insight into how both 
phase transitions are connected. 

After a short introduction we will discuss several aspects of this new relation between the 
Polyakov loop and the eigenvalues, in particular we focus on the volume scaling and the continuum 
limit. We will present numerical results for both quenched and dynamical QCD and will also 
compare to the free case. 



2. Introduction 



Starting point of our discussion are the eigenvalues A of the massless staggered Dirac operator, 
D*F = ibiA*F with Xreal > 0. The massless staggered Dirac operator is defined by 



xy ■ 



1 ^ 



(2.1) 



where r\xii is the usual staggered phase factor and Ux^i are the link variables. Note that we use 
periodic boundary conditions in all four directions for the calculation of the eigenvalues. The 



Polyakov loop is defined by P = ^^-^ Y.n Tr^ 
the eigenvalues in the following way 



and can be expressed in terms of 



frJlZ y [\.X^' +7-X^' +7* 



A, 



(2.2) 



Ns and Nt is the spatial and temporal extension, respectively, and Z3 = {l,z,z*}. Note that Nt has 



to be even for staggered fermions. The sum over / in Eq. ^^is meant to sum over all eigenvalues 
where stands for the eigenvalues calculated on a given gauge configuration which is Z3- 
rotated by X € Z3. To be less confusing, for a given gauge configuration we generate all three 
Z3-rotated gauge configurations and calculate all eigenvalues of the Dirac operator for each of 
the three configurations. The Polyakov loop can then be expressed as a sum over all eigenvalues 
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Figure 1: On the left hand side we plotted P(A) in the complex plane for a typical gauge configuration 
at r > Tc for A', = 4. The right hand side shows a sketch of the corresponding eigenvalues in the real and 
complex sector of P. 



calculated on all three rotated gauge configurations, Eq. g^. By looking at Eq. one immediately 
may ask what part of the eigenvalue spectrum contributes most to the Polyakov loop. To answer 
this question one introduces the following cumulative sum, 

22^' ^ 



l-A,^'+z-A,/^. 



(2.3) 



where we sum over all eigenvalues up to a certain (maximal) value A. Let us briefly make some 
comments on this formula before discussing an example. Note that the sector where the Polyakov 
loop sits for a given configuration (at T > Tc) is solely determined by multiplying the with the 
appropriate Z3-factors. Let us assume for the moment that the gauge configuration corresponding 
to has P in the real sector. Remember that the two complex sectors of P are physically equiv- 
alent. Therefore, the eigenvalues in the complex sectors of the Polyakov loop are approximately 
the same. A,- j « A, Making use of this we obtain for Eq. 2.2 the following approximate relation, 
P{^) ~ Ia,,x<a(^,>o/ - ^complex)- TWs relation illustrates the fact that P(A) is built up by the 



"response" of the eigenvalues on the different Z3 sectors. As an example we have plotted P{)i) for 
a typical gauge configuration in the complex plane on the left hand side of figure [T[ On the right 
hand side we plotted a sketch of the corresponding eigenvalues calculated in the real and complex 
sector of the Polyakov loop. We observe that for small eigenvalues the complex sector dominates. 
By looking at our previously derived relation we see that this results in negative values of P{X) 
while for large eigenvalues it is the other way round. This shows that the change of the eigenvalues 
with respect to the different sectors of the Polyakov loop is crucial. 

Finally, we perform the ensemble average on the absolute value of P(A), 



(2.4) 



This is the object we will study for the rest of our discussion. At this point, let us draw the reader's 
attention to Refs. [Q-^ were similar investigations has been performed using staggered and Wilson 
fermions. 
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Figure 2: Results for (|P(A)|) for our quenched configurations. The upper (lower) plots correspond to 
T > Tc (T < Tc). The plots on the left hand side illustrate the volume dependence, the plots on the right hand 
side are rescaled such that {\P{A,max)\) is the same for all three volumes. In the later case we observe that for 
each j3 all three curves lie above each other telling us that ( ) | ) has the same volume scaling than P. 

3. Data 

In this section we present several numerical results for {\P{X)\). We start with results for the 
quenched case where we have used standard Wilson gauge action. The statistic varies from ~ 10 
configurations for the largest lattice 8*^ up to ~ 100 configurations for the smallest lattice. The 
eigenvalues were calculated on a single work station using the ARPACK library |^. 

Let us first take a look at the plot on the upper left side of figure ^ where T > Tc and the 
Polyakov loop is finite. What we notice right away is that the main contribution to the Polyakov 
loop comes from the large eigenvalues which is somewhat surprising since the physically relevant 
part of the spectrum should be the infrared. We will comment on this later on. Another surprising 
observation is the dip of the curves at A 1.7 — 1.8. Naively, one might expect that the cumulative 
sum {\P{X)\) is a monotonically increasing function. However, by looking at figure |I| this behavior 
becomes clear. The dip in (|P(A)|) (curves in the upper left plot of figure ^ corresponds to the 
region where (see figure [l]) passes zero, the bump in the curves at A 1.4 corresponds to the 
region where P{X) takes its negative values. This structure seems to be quite interesting and one 
may ask whether it will survive the infinite volume and the continuum limit. 
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Figure 3: The upper plots show the behavior of (|f (A)|) as we decrease the lattice spacing a at fixed 
T « l.2Tc. The plot on the right hand side is plotted in physical units. The lower plot shows ( |P(A ) | ) in the 
free case for the same values of Nt at quite large spatial volume. 



In figure ^ we illustrate the volume dependence of (|P(A)|). On tlie upper plots we sliow 
results for three volumes 4^, 6^, 8^* with Nt =4 at a. temperature slightly above Tc. The plot on the 
upper right hand side shows rescaled curves where we have fixed {\P{X)\) = (1^83x41) at A = Xmax 
for all three volumes. Remember that P{X,„ax) is just the ordinary Polyakov loop P. Beside the 
curve corresponding to the smallest volume, which shows small deviations, the curves lie above 
each other. This observation tells one that for T > (for sufficiently large volumes) (|P(A)|) has 
the same volume scaling than the Polyakov loop itself. This in turn means that the structure will 
survive the infinite volume limit. 

For T <Tc the situation is similar. On the lower left hand side of figure |we plotted (|P(A)|) 
for the same three volumes. The lower right plot shows the corresponding rescaled curves where we 
again fixed (|P(A,„a,v)|) = (|-f83x4l)- Again, the curves lie above each other showing that (|P(A)|) 
scales like the Polyakov loop also below T^. Because below the Polyakov loop vanishes in the 
infinite volume limit we find that also (|P(A)|) will vanish in this limit. So in this case the structure 
does not survive the infinite volume limit. Note that knowing the volume dependence of (|P(A)|) 
will keep the computational costs significantly lower because one do not have to perform expensive 
computations on large volume lattices. 
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Figure 4: We plotted (|P(A)|) using dynamical (j3 = 3.26,3.36) and quenched (j3 = 5.4,5.8) configurations. 
The dynamical configurations were generated with masses m^=0.0065 and mi=0.065. The temperature in 
the dynamical case differs slightly from the quenched one. 



Let us now look at how (|P(A)|) behaves as the lattice spacing a — > 0. Since for T < T^- 
vanishes in the infinite volume limit anyway we will discuss only results for T > Tc. The 
upper plots in figure |3| illustrate our results for Nt =4,6,8 at fixed T ^ \ .2Tc plotted against A in 
lattice and physical units. Let us focus on the upper left plot. We find that for Nt = 6 there is a 
bump and a dip in the curve similar to that in the curve for Nt = 4. For N = S we also observe 
a large bump but because in this case the statistic is quite limited and the spatial volume is rather 
small the signal is quite noisy. We notice that, as we go to smaller lattice spacing, the structures 
in the curves move towards the ultraviolet (UV) part of the eigenvalue spectrum and the Polyakov 
loop obtains its final value more from the very end of the UV part of the spectrum. 

Let us compare these results to the free case where an analytical expression for the eigenvalues 
is known. The lower plot in figure || shows our results in the free case for the same three values of 
Nt at quite large spatial volume. We find that the shape as well as the position of the bumps and 
dips of the curves are surprisingly similar to the corresponding results of the quenched data (upper 
left plot). This leads to the following conclusion. As we approach the continuum limit at 
small values of A is essentially zero. At somewhat large eigenvalues (|P(A)|) starts to show wild 
fluctuations which cancel out at the very end of the UV part of the eigenvalue spectrum where the 
Polyakov loop obtains its final value. 

We remark that it might be not too surprising that the Polyakov loop is dominated by the UV 
part of the eigenvalue spectrum as the Polyakov loop is related to the propagation of an infinitely 
heavy quark. By looking at the quark propagator in the spectral representation, 

S(x,y)=l^ — 7—. , (3.1) 

where i//;^ i^) the normalized eigenvectors of the Dirac operator, we note that the eigenmodes 
under the sum are weighted by (A + im)^K For a very heavy quark (as m — oo) the relative weight 
of each eigenmode becomes approximately the same. Therefore, UV eigenmodes can dominate the 
propagation of an infinitely heavy quark, i.e. the Polyakov loop. 
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Finally, we compare the quenched results to dynamical results using p4fat3 fermions with 
improved gauge action and quark masses m^=0.0065 and =0.065, see Ref. [0]. On the left (right) 
hand side of figure^ we plotted (|P(A)|) for T > Tc {T < T^). Surprisingly, there is no qualitative 
difference in the behavior of (|P(A)|) in the dynamical case. 



4. Summary 

In this work we have studied the connection between the Polyakov loop and the eigenvalues of 
the Dirac operator using {\P{X)\). We have focused on the volume dependence and the continuum 
limit. We have found that the dominant contribution to the Polyakov loop comes from the very 
end of the UV part of the eigenvalue spectrum. We also compared our results to the free case. A 
comparison between full QCD and quenched QCD seems to show no qualitative difference. 

Our findings suggest that the dependence of the eigenvalues on the different Z3 sectors of the 
Polyakov loop seems to be crucial. Aiming at the connection between confinement and chiral sym- 
metry breaking our findings can be concluded in the following picture. Above Tc the eigenvalues 
A show a strong dependence on the different sectors of the Polyakov loop resulting in wild fluctu- 
ations in These fluctuations cancel out in a way that the Polyakov loop obtains its finite 
value from the very end of the UV spectrum. At the same time, since the chiral condensate is zero 
above T^ the infrared (IR) part of the spectrum shows a vanishing density of eigenvalues. 

Below Tc the dependence of the eigenvalues on the Polyakov loop sectors vanishes in the 
infinite volume limit which leads to a vanishing (|P(A)|). In particular (|P(A)|) vanishes at the UV 
leading to a Polyakov loop (|P|) =0. At the same time, the eigenvalue density at the IR part of the 
spectrum becomes finite since chiral symmetry is broken. 
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